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Abstract 

The concept of ratchets, driven asymmetric periodic structures giving rise 
to directed particle flow, has recently been generalized to a quantum ratchet 
mechanism for spin currents mediated through spin-orbit interaction. Here 
we consider such systems in the coherent mesoscopic regime and generalize 
the proposal of a minimal spin ratchet model based on a non-interacting clean 
quantum wire with two transverse channels by including disorder and by self- 
consistently treating the charge redistribution in the nonlinear (adiabatic) 
ac-driving regime. Our Keldysh-Green function based quantum transport 
simulations show that the spin ratchet mechanism is robust and prevails 
for disordered, though non- diffusive, mesoscopic structures. Extending the 
two-channel to the multi-channel case does not increase the net ratchet spin 
current efficiency but, remarkably, yields a dc spin transmission increasing 
linearly with channel number. 

Keywords: ratchets, spin electronics, mesoscopic quantum transport 
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1. Introduction 

The appealing physical concept of converting energy from randomly mov- 
ing Brownian particles in asymmetric set-ups into a directed particle flow, 
possibly against an external load, has led to an enormous amount of works 
establishing an own field at the interface of transport and nonequilibrium 
statistical physics. Peter Hanggi, as one of the founders of this branch of 
statistical physics, coined the term "Brownian motors" for such systems 
Ratchets, spatially periodic structures with broken left-right symmetry op- 
erating far from equilibrium and thereby generating directed particle motion 
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in the presence of unbiased time-periodic driving constitute one important 
class of such systems. The ratchet mechanism was first discovered for classical 
Brownian particles j2|, [s], 0] and then generahzed to quantum dissipative sys- 
tems [5|] . Later this concept has been extended to the coherent regime where 
corresponding ratchets and rectifiers have gained increasing attention, in par- 
ticular after the experimental demonstration of ratchet-induced charge flow 
in periodically arranged lateral quantum dots based on a two-dimensional 
semiconductor heterostructure p]. Coherent rectifiers are characterized by 
phase-coherent quantum dynamics in the central periodic system in between 
leads where dissipation takes place. For a recent comprehensive review of 
the whole, broad field of ratchets and artificial Brownian motors see 0]. 

While nearly all the works in this field have addressed the problem to 
achieve unbiased directed particle transport, we have generalized this con- 
cept to the notion of spin ratchets, corresponding set-ups which allow for 
generating spin currents, partly even in the absence of charge currents. In 
this respect, "Zeeman ratchets" that are based on an asymmetric, spatially 
periodic magnetic field are closest to the usual charge ratchets: Owing to 
the Zeeman term in the Hamiltonian, spin-up and -down electrons experi- 
ence opposite asymmetric periodic potentials, which should give rise to a net 
flow of charge carriers with different spin polarization in opposite directions, 
corresponding to a pure spin current j8| . Similar results have been found for 
a one-dimensional quantum wire with strong repulsive electron interactions 
[oj. However, contrary to particle ratchets with preserved particle number, 
spin-p olarizat ion is a volatile property and subject to spin relaxation. Hence 
in 10| we showed, both conceptually and numerically for a realistic model 
of a Zeeman ratchet including spin-flip processes that pure spin current gen- 
eration in coherent mesoscopic conductors is indeed possible. In parallel we 
devised the concept of a "spin-orbit ratchet" [IH, a setting where the spin 
orbit interaction (SOI) in a quantum wire is employed to generate a spin cur- 
rent. Contrary to particle ratchets, which rely on asymmetries in either the 
spatially periodic modulation or the time-periodic driving, a spin-orbit (SO) 
based ratchet works even for symmetric electrostatic periodic potentials, due 
to the spin-inversion asymmetry of the SOI. As a result, no ratchet charge 



current is produced in parallel, leading to a pure ratchet spin current [11 
As possible realizations we have in mind systems based on semiconductor 
heterostructures with Rashba SOI 12| that can be tuned in strength by an 
external gate voltage allowing to control the spin evolution. 

Such a SO-based ratchet is an example of a device for "mesmerizing" 
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semiconductors [13[, i.e. for magnetizing semiconductors without using mag- 
nets, in contrast to the usual method of spin injection through ferromagnets. 
Also in this sense spin ratchets have much in common with spin pumping 
i. e. the generation of spin-polarized currents at zero bias via cyclic parameter 



variation. Different theoretical proposals based on SO and Zeeman [15 
mediated adiabatic spin pumping in non-magnetic semiconductors have been 
suggested and, in the latter case, experimentally observed in mesoscopic cav- 
ities [igI]. Adiabatically driven spin ratchets, however, also differ from adia- 
batic spin pumps, since ratchets operate with a single ac driving parameter 
in the nonlinear bias regime, while the usual proposals for adiabatic pumping 
include a cyclic variation of at least two parameters at zero external bias [17 . 

Recently, the concept of SO-based spin ratchets has been generalized to 
the quantum dissipative regime. To this end the above setting was extended 
by coupling the orbital degrees of freedom additionally to an external bath 
(within a Caldeira-Leggett model H) representing, e.g. effects of phonons. 
While SO-mediated spin-phonon coupling usually leads to spin relaxation, it 
could be shown that for this ratchet set-up the opposite is true: a finite, pure 



spin current is generated [22|, |23| . This is remarkable as it means that ther- 
mal energy from the bath is converted, via the SO coupling, into a directed 
spin current with aligned charge carrier spins; in this sense such a dissipa- 
tive spin ratchet can be viewed as a "Brownian spin motor". An extension 
to an additional in-plane magnetic field which allows for further controlling 
the spin current can be found in [2J|. The study of dissipative SO ratchets 
took another twist when it turned out that a finite charge current arises in a 
SO ratchet even if both the spatially periodic potential and the time-periodic 



driving are symmetric [25| . While directed transport is well known to appear 



for symmetric potentials as long as the driving contains a time asymmetry 



26|, |27|, |28| , it seems paradoxical that one can even go without this symmetry 
breaking. The charge ratchet mechanism for this space- and time-symmetric 
case results from the interplay between quantum dissipation and SO-induced 
spin-flip processes and thereby from a hidden symmetry breaking through the 
SO coupling. This new class of space- and time-symmetric charge ratchets is 
hence built on the intrinsic spin nature of the particles to be transported, ex- 
hibiting certain conceptual similarities with the classical "intrinsic ratchets" 
discussed in [291]. 

In the present paper we however focus on coherent SO ratchets and extend 
the minimum model introduced in [ll| in several directions in order to gain 
an improved and systematic understanding of the spin ratchet mechanism 
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and its limits for set-ups describing realistic mesoscopic devices. 

First, most investigations of ratchets disregard interaction effects and 
usually model the time-periodic driving in terms of the bare homogeneous 
external field. However, in particular for charge ratchets in the nonlinear 
bias regime, the actual voltage drop across the system may deviate from 
simple bias models often used, e.g. from a linear model, and hence can affect 
the ratchet current which is known to depend sensitively on the details of 
the system. Here we will consider in a self-consistent treatment the charge 
redistribution in the ratchet due to the external bias giving rise to a nonlinear 
voltage drop and thereby presumably altering the resulting spin-dependent 
transmission through the device. 

Our study includes, second, an analysis of disorder effects usually un- 
avoidable in mesoscopic devices. It is well known that impurity scattering 



in a SO medium usually leads to Dyakonov-Perel-type spin relaxation |30 
This mechanism will obviously counteract the generation of spin-polarized 
currents, and hence it is important to understand to which extend the spin 
ratchet mechanism is robust against disorder effects. 

Third, as shown in at least two SO-coupled transverse modes in a 
quasi-one- dimensional wire with at least one electrostatic potential barrier 
are required for generating a net spin current. While an increase in the bar- 



rier number on the whole leads to an enhanced ratchet spin current [11 
the spin current dependence on the number n of transverse modes, respec- 
tively the wire width, remains to be investigated. Our analysis shows that, 
on average, the ratchet spin current does only marginally increase with n. 
However, interestingly, we find for a given (dc) bias a linear increase of the 
spin-polarized transmission with n. 

The paper is organized as follows: After introducing the system and 
the numerical method in Sec. 2 below, in Sec. 3 we will first consider spin- 
polarization effects in dc-transport through two-dimensional ribbons with 
Rashba SOI and outline an important underlying polarization mechanism. 
In Sec. 4 we then study the SO ratchet response for the nonequilibrium and 
disordered case and conclude with a number of remarks in Sec. 5. 



2. Outline of the system 

We consider a quantum wire (oriented in x-direction) which can be re- 
garded as being realized in a two-dimensional electron gas (2DEG) in the 
(x, |/)-plane. A typical structure is visualized in Fig. [1^. In the central region 
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of the wire Rashba SOI is present, which is described by the Hamiltonian 



Hr = ^a{x)a^py - ^(Ty [a{x)px + Px(y{x)] . (1) 

To avoid unwanted reflections at the interfaces between the SO-free leads 
and the scattering region with finite SOI a{x), we adiabatically turn on the 
Rashba SOI (see Fig. Wp)- Since Rashba SOI is only present in the central 
region of the wire [sij], we do not face difficulties with the definition of a spin 
current inside leads with SOI (32| . 

In addition we consider A^b identical potential barriers, which are located 
in the wire (see Fig. They are modeled by the electrostatic potential 



f4arr(a;) 



if/B |l - COS [l^ {x + INbLb)^ } for \x\ < NbLb/2 

elsewhere 



(2) 



The Hamiltonian of the whole system then reads 

2 2 

H = + Uconiix, y) + [/barr(x) + 6U,,{x, y) + Hr , (3) 

2m* 

where Uconi{x,y) is the hard-wall confinement potential and 6Ues{x,y) is an 
electrostatic potential due to the rearrangement of charges at finite bias, see 
Sec. 4.1. 

A quantitative study of nonequilibrium and disorder effects on ratchet 
transport requires a numerical approach. For the respective transport calcu- 
lations, we use a tight-binding version of Eq. ([3]), i.e. the system is discretized 
on a square grid with lattice spacing a. By using an efficient recursive al- 
gorithm for computing the lattice Green functions [s^ we determine the 
relevant spin-dependent transport properties of the system. From the nu- 
merically obtained S'-matrix elements we then calculate the spin resolved 
quantum transmission probabilities T^^, between the left and right lead. 
Here T^fi describes the probability for an electron with spin state a' to be 
transmitted from the left entrance lead into the right exit lead with spin 
state a. The spin state a = ± is defined with respect to a spin quantization 
axis pointing in ^/-direction. Such an in-plane polarization is often referred 



to as Edelstein effect [3J] • From the spin-resolved transmissions we can then 
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calculate the total and spin transmission, respectively: 

^(^) = (4) 

(T a' 

TsiE) = $]K^,(i5)-T«^,(i?)] . (5) 

a' 

Contrary to the total transmission, Eq. which in a two-terminal setup 
is symmetric with respect to interchanging the leads, i.e. T = T^^ = T^^, 
the spin transmission (and naturally also the associated spin current) can 
differ from lead to lead: in general T^^ ^ T^^ owing to the spin inversion 
asymmetry of the SOI. Therefore it is necessary to specify the lead where the 
spin current is evaluated. We will calculate the spin transmission from the 
left to the right lead using the abbreviations Ts{E) = T^^{E), see Eq. (|5]), 
and Taai{E) = T^^{E) employed below. For later use we furthermore define 
dimensionless energies and Rashba SOI strengths (denoted by a bar): E = 
[h'^/{2m*a'^)]E and a = {m*a/h^)a. 
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Figure 1: a) Spin ratchet setup consisting of a quantum wire along the x-direction. The 
width of the leads is W while the central region has a width Wq. b) strength of the Rashba 
spin-orbit interaction a(x); c) electrostatic potential barriers J7barr(a;); d) additional elec- 
trostatic potential offset Uos{x). 
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3. Spin polarized dc-transport 

Before turning our attention to tlie ratcliet behavior upon ac-driving we 
investigate the dc-transport properties of the system outlined in Fig. [T] in 
order to demonstrate the effect of the SOI underlying the spin ratchet mech- 
anism. To this end we evaluate the charge/spin currents in the leads in 
response to a fixed finite, applied bias making use of the expressions from 
the Landauer-Biittiker formalism. The charge current for coherent transport 
in a quantum wire reads 



Here f{E; /i) is the Fermi-Dirac distribution function and /xl/r is the chemical 
potential of the left / right lead. Correspondingly, the spin current in the right 
lead reads 



In Fig. [2] we show the numerically computed transmission probabilities T{E) 
and Ts{E) in linear response for a system with five potential barriers. Fig. [1] 
As shown T{E) (full black line) exhibits two combs of four transmission 
peaks each due to resonant tunneling through the array of potential barriers, 
reflecting precursors of minibands arising for an infinite array of barriers. 
The two transmission combs belong to transport of states with transverse 
mode number n = 1 and 2, respectively. Furthermore, the system exhibits 
strong spin polarization in +y direction, i.e. Ts > (dashed red line), for a 
wide range of Fermi energies above the subband energy of the second channel. 
This indicates that SO mixing of at least two modes is required for a nonzero 
spin transmission. 

We now outline the basic mechanism which causes such a spin polariza- 
tion. To this end we employ a Landau-Zener model for a single barrier which 
was first introduced by Eto et al. [sl] to describe spin-dependent transport 
across a quantum point contact (see also [s^, [s^]). In Fig. [3] we show the SO- 
split parabolic dispersions (in fcaj-direction) of charge carriers in the two first 
transversal subbands of the quantum wire, relative to the fixed Fermi energy 
Ep (horizontal dashed line), at three different positions A-C of the potential 
barrier (along the x-direction). The SOI further couples different parabola 
branches and leads to a small avoided crossing (position marked by a dashed 
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Figure 2: Charge (full black line) and spin (dashed red line) transmission probability in 
linear response for the system depicted in Fig.[TJ Parameters used are: Lb — 10a, Nb = 5, 
ao = 0.15, W = 25a, Wo = 15a, Ub = 0.2, Uos = 0.15. 



box) between the states (n, cr) = (2, +) and (1, — ) (full red and dashed blue 
line) as shown in Fig. [31 In [35|] it was found that upon traversing the barrier 
(A— J-B— j-C), first higher transversal modes of the wire become depleted (see 
position B in Fig. [3]). After passing the barrier top the SOI gives rise to a 
spin- dependent repopulation of these higher modes between position B and C 
when the Fermi energy passes adiabatically the afore mentioned anticrossing. 
For the simplest situation of two occupied transversal modes and a single 



barrier (shown in Fig. [3]) the spin transmission can be estimated as 36 

Tg = 2P , where P = 1 - e'^ (8) 

is the transition probability between subbands with different spin polar- 
ization {n = l,cr = — ) ^ {n = 2,a = +)• This quantity was evaluated in 



Ref. [35| using Landau-Zener theory [38|, |39|, where A > parametrizes the 



adiabaticity of the transition. In the diabatic limit, A — )■ 0, Landau-Zener 
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Figure 3: Sketch of the parabolic energy dispersion for two transverse modes of the 
quantum wire in the presence of an electrostatic barrier shown at three different posi- 
tions (A)-(C). The spin-splitting due to SOI is indicated by the shift of the parabolas 
in ±A:2:-directions. The dotted square indicates an avoided crossing of transversal modes 
(n, a) = (2, -|-) and (1, — ) with different spin polarization, where the spin-flips can happen, 
see Eq. ([8)). The Fermi energy is marked by the horizontal dashed line. 

transitions preserving the state (1, — ) dominate, while in the adiabatic hmit 
the (1, — ) state (dashed blue line in Fig. [3]) changes its character into (2, +) 
(solid red line). The Landau-Zener parameter A depends on the form of 
the barrier, the confinement potential and the SOI. As two of us showed in 
Ref. ^4^, Ts increases with the length of a point contact constriction, which 
corresponds to the adiabtic limit of an increasing length Lb of the barrier in 
the set-up investigated here. 

Generalization of these results to a higher number n of contributing transver- 
sal channels shows that this spin polarization effect is not limited to the two- 
channel case. In Fig. IHwe present our numerical results for T and Ts for fixed 
Fermi energy as a function of the width W of the conducting stripe which is 
proportional to the number of occupied transversal modes: n = kplV/ir with 
Fermi wave number kp. We see that both the total and the spin transmission 
increase linearly with W yielding a constant polarization ratio T^/T of the 
transmitted electrons. This ratio can be controlled by the strength of the 
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n 



Figure 4: a) Total transmission (full black line) and spin transmission (dashed red) in 
linear response as a function of the channel number n in the lead for a scaled SOI strength 
cto = 0.15. b) Corresponding ratchet spin transmission ATs, Eq. (fT4|) . for a linear voltage 
drop at finite bias Uq — 0.02. Inset of panel a): Spin polarization ratio Ts/T as a function 
of the scaled SOI strength at fixed width W = 250a corresponding to n = 56 transverse 
channels. Common parameters: Lb = 60a, A^b = 1, E-p ~ 0.36, Ub — 0.22, Uos ~ 0. 



SOI as shown in the inset of Fig. HI 

While the hnear increase of T for large W is expected and can be straightfor- 
wardly understood from the increasing number of open channels carried by 
the leads, the linear growth of Ts with wire width is remarkable. In order to 
explain this linear rise of Tg, we have to consider the number of SO-coupled 
channels participating in the spin polarization mechanism and to general- 
ize the polarization mechanism described above to the case of sequences of 



Landau-Zener-type transitions (see j4l| for an explicit treatment of the case 



n = 3). Without invoking a further detailed analysis of the SO-induced 
coupling mechanisms, a simple density-of-states argument shows that the 
number of relevant states in a critical energy window of order of the barrier 
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height below the fixed Fermi energy scales also linearly with W. We finally 
note that, while the l^-linear behavior of T is also obtained classically, the 
linear increase of Ts holds true for fully coherent, ballistic transport and pre- 
sumably does not prevail for widths larger than the phase coherence length. 
We have not explored its dependence on disorder. 

4. Spin ratchet effect 

After having explained the main mechanism responsible for spin polar- 
ization in dc-transport across a single barrier, we are now prepared to turn 
our attention to the more general case of a periodic arrangement of barri- 
ers and ac-driving, i.e. the study of the spin ratchet effect. We consider an 
adiabatic ac-driving assuming that the external bias /xl — A^r is varied on 
a timescale much longer than the relevant timescales for electron transport 
through the system, which is basically the dwell time of the electrons in the 
ratchet scattering region. In experiments, the efficient operation of coherent 
charge ratchets in this adiabatic regime has already been confirmed j^. 

To be specific, here we consider an adiabatic square- wave driving with 
period to, where the chemical potential of the left /right reservoir. 



is periodically switched [fiL/R(t) = fJ'L/nit + to)] between the two rocking 
situations with bias difference ±Uo = /xl — /^r with Uq > 0. For this adiabatic 
square wave driving, the ratchet is in a steady state in-between the switching 
events. Thus, we can use the respective expressions and ([7]) for the dc- 
charge and spin current to calculate the averaged spin ratchet currents. For 
the driving of Eq. iQ the ratchet currents are given by the average between 
the two rocking conditions +Uo and —Uq: 




for < t < to/2 , 
for to/2 <t<to, 



(9) 




(10) 



{I^{Ep,Uo)) = ^[I^{Ep,+Uo) + IHEf,-Uo)] 



(11) 
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where 



A/(E;Ep,f/o) = f{E;Ep + Uo/2)-f{E;Ep-Uo/2), (12) 
ATiE;Uo) = TiE;+Uo)-TiE;-Uo), (13) 
ATs(E;f/o) = Ts{E;+Uo)-TsiE;-Uo). (14) 

In linear response [6Ues{x,y) = in Eq. ([3])] the quantities AT{E; Uq = 0) = 
ATs{E;Uo = 0) vanish and thereby also the charge currents, Eq. (ITUI) and 
spin ratchet currents, Eq. f llip . This implies that the system has to be driven 
into the nonlinear regime, typical for ratchets. Therefore, from now on we 
apply a finite driving bias Uo to operate the spin ratchet. This bias deter- 
mines the voltage drop 6Ues{x, y) in Eq. (|3]), which describes the electrostatic 
potential due to the rearrangement of charges in the wire compared to the 
linear response case. 

The spin ratchet mechanism, resulting from a finite ATs{E; Uq), can be 
qualitatively understood from the Landau-Zener model used in Sec. 3 to 



explain the spin polarization mechanism of a single barrier. In [ll[, the 
expression (|8]) for the spin-flip probability was extended to include a finite 
voltage drop: 

^ ^ ^ " ''''^ { {d/dx){U^Jix) + 6UUx)) } ' ^^^^ 

where rj depends on the SOI strength and the confinement potential of the 
quantum wire, and the derivative is evaluated at the position x of the avoided 
crossing marked by the dashed box in Fig. [31 We see that P depends on the 
amplitude of the driving voltage via the gradient (9/9x)([/barr(a;) + SUcs{x)) 
yielding different transition probabilities for forward and backward bias. 
Since the transitions are induced when the degeneracy points (marked by 
the square in Fig. [3]) cross the Fermi energy, the value of d{6Ues)/dx in 
the vicinity of the barriers is important for the appearance of the spin 
ratchet effect. Therefore, this model predicts finite ATs for finite Uq, since 
P{+Uo) 7^ P{—Uq). In summary, the spin ratchet effect predominantly re- 
sults from the deformation of the potential barrier due to voltage drop SUes 
which enters into the SO-mediated spin-flip processes. 

We flrst examine the ratchet spin current dependence on the number of 
transversal modes by approximating 6Ues{x, y) by a linear function. In Fig.Hh 
we see that, in analogy to the dc-transport quantities T and Ts (shown in 
Fig. H^), the ratchet spin transmission ATs also exhibits on average a linear 
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increase for large n, though with a very small slope. Its overall magnitude 
stays below one, and it exhibits larger values at small n due to the strong 
fluctuations in that regime. We can conclude that increasing the number of 
transverse modes does not significantly enhance the ratchet spin current. 

4-1 ■ Self- consistent calculation of the voltage drop 

In most models for ratchets, the effective voltage drop 6Uesix,y) along 
the wire has been approximated by a linear function. In this section we will 
investigate the role of the voltage drop by comparing spin ratchet signals 
resulting from 6Ues{x,y) with its linear approximation. To obtain 6Ucs{x,y) 
in the scattering region, we self-consistently solve the Schrodinger equation 
and the Poisson equation. To this end we adopt the approach introduced for 
ratchets in fS], which we describe in the following. Within this approach 
we absorb all electrostatic potentials (e.g. the potential due to donor atoms), 
which do not change upon variation of system parameters, as do the Fermi 
energy or the bias voltage, in the confinement potential of the wire. Then 
we only have to consider the rearrangement of the electrons, 6n = n — hq, 
due to a finite bias in order to determine the voltage drop in the system. 
For the calculation of the electron density we determine the lesser Green's 
function of the system via the Keldysh equation 43 



(16) 



where Q'' is the retarded Green's function of the system, which we evaluate 
via a recursive Green's function method 4^. Furthermore, the lesser self- 
energy can be expressed in terms of the retarded self-energies of the 
individual leads i, 



i 

Finally, the electron density of the system is given by 

i ("^ 

-— / dETi[g<{r,r-E)]. 



n r 



We then obtain the change in the electrostatic potential, 51]^, 
by solving the corresponding Poisson equation 44 



n[r) — nQ[r) 



(17) 
(18) 

(19) 
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where is the vacuum permittivity and Er is the material specific relative 
static permittivity. For the evaluation of 6Ues{r) it is useful to distinguish 
the contributions from the leads, 6Uiea,d{r), and from the rearrangement of 



the electrons in the scattering region, 5Usr{r) |44j . |45 



SU^f^ = 5?7iead(r) + SU^f^ (20) 

Here, SU\cad{r) solves the Laplace equation V^5[/iead(^ = with the bound- 
ary conditions SUiead{x = ±L/2) = and is therefore given by a linear 
function between both contacts [i^. On the other hand, 5U^^{r) solves the 
Poisson Eq. (fT9|) with boundary conditions 5Usv{x = ±L/2) =0. To obtain 
SUsrir) from this equation we assume that the electron density inside the 
leads is much higher than in the scattering region. In the simulations we 
realize this by introducing an additional electrostatic offset potential in the 
scattering region, see e.g. Fig. [T]i). As a result the electrostatic potential 
profile close to the leads is fiat, i.e. n(f^ fa ^^-0(^5 which enables us to calculate 
SUsr from the Poisson Eq. f|T9|l with vanishing SUsr for |r| — )■ 00, yielding 

6UAf) = ^ [ dV'^^^^^l^^ . (21) 
AnErSo J \f-r'\ 

Now we can compute the electrostatic potential SUes{x,y) at finite bias volt- 
ages. To this end, we start with an initial guess for 6Ues{x,y) and calculate 
the electron density n{x, y) for this case via Eq. ffTSjl . By solving Eq. fl2Tl) we 
obtain a new potential 6Ues{x, y), which in turn can be used to calculate the 
corresponding density n{x,y). This procedure is repeated until convergence 
is reached. In practice, we do not directly iterate between Eqs. flTS]) and fl?I]) . 
but we use the so-called Newton-Raphson method, which has been success- 
fully applied to similar non-equilibrium problems [46[ 47 1 and significantly 



improves the convergence of the self-consistent calculations. 

Before turning to the spin current calculation we first present a symmetry 



analysis of the spin- resolved transmission probabilities [10|, |48[ . This is help- 
ful to simplify the expressions for the ratchet currents. To be specific, the 
Hamiltonian ([3]) is invariant under the operation of 

V = -iCRuR^a, , (22) 

since the symmetry relations 

6Ucs{x) = -5f/cs(-a;), f/barr(a;) = t/'barr(-a;), a{x) = a{-x) (23) 
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are fulfilled. In Eq. ( 12^ . inverts the x-coordinate, switches the sign of 
the bias voltage (±?7o ^ TUq) and C is the operator of complex conjugation. 
The first equality, 6Ues{x) = —6Ues{—x), is due to the spatial symmetries of 
the system and the use of the same computational scheme for both forward 
and backward bias. In Qdl we showed that due to this invariance the relation 



(24) 



T,y{E,±Uo)=T,'AE,TUo) 



is fulfilled. As a consequence, the charge ratchet current (J), Eq. f[TU|) . for 
this system vanishes. On the other hand, the expression for the spin ratchet 
current, Eq. (fTT]) . can be simplified through the above symmetry relation: 



(I'm) 



1 



dEAf{E;Uo) T+,_(E,+[/o)-r_,+ (^,+f/o 



(25) 



Hence, it is sufficient to calculate the spin-fiip transmission probabilities T_^- 
and T_ + for a single rocking condition. 

In Fig. Owe present the calculated spin ratchet conductance (I^) /Uq as a 
function of the amplitude of the driving voltage Uq, see Eq. ([9]). As the 
2DEG material we choose InAs with the parameters m* = 0.024 mg, g* = 15 
and Ej- = 15.15. Furthermore, the lattice spacing is set to a = 10 nm. We 
find a finite spin current with a direction depending on the Fermi energy. We 
compare the results of a linear voltage drop between x = —L/2 and x = L/2 
(dashed curves in Fig. |5]) with those obtained from the self-consistent calcu- 
lation of 6Ues{x,y) (solid curves) for several representative Fermi energies. 
The self-consistent calculations and the linear voltage drop model yield sim- 
ilar results for (/^) except for the case Ep = 6.03 meV, where transport is 
dominated by resonant tunneling and hence depends sensitively on the de- 
tails of the potential. In Fig. [H] we depict the spatial distribution of the self- 
consistently calculated voltage drop 6Ucs along the ratchet wire. We see that 
the linear ramp is indeed a good approximation for the three higher Fermi en- 
ergy values considered. This also explains the good qualitative agreement in 
Fig. [5] between the spin conductance for the self-consistently determined 6Ues 
and the linear voltage drop model, respectively. For the case of Ep = 6.03 
meV the resulting voltage drop shows the most pronounced non-monotonic 
behavior in Fig. |6l This causes a misalignment of the energy levels in the po- 
tential valleys and therefore a reduction of the miniband-mediated resonant 
transport. As a consequence the spin ratchet current is overestimated in the 
linear voltage drop model for this value of Ep. 
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Figure 5: Ratchet spin conductance {I^)/Uo as a function of the driving voltage Uq for 
four different values of Ep. We compare the results for a linear voltage drop (dashed 
lines) with those for the self-consistently determined SUes{x,y) (solid lines). Parameters: 
Lb = 100 nm, Nb = 5, ao = 4.76 • 10-"eVm, W = 250 nm, Wo = 150 nm, Ub = 3.17 
meV, Uos = 2.38 meV. 



4.2. Influence of disorder 

So far we have studied the case of disorder free, clean conductors. How- 
ever, in reahstic experimental samples dopands, crystal defects and impurities 



give rise to momentum scattering. This in turn can cause spin relaxation [30 



which might limit the performance of the spin ratchet. In order to investigate 
the role of impurity scattering on the ratchet effect we again consider the de- 
vice shown in Fig. [1] but with W = Wq = 15a and without the additional 
potential offset shown in Fig. [T]i. For the sake of computational feasibility in 
the following we assume a linear voltage drop in the central region, since we 
have seen that such a model represents a fair approximation for the actual 
voltage drop for a wide parameter range, and for the sake of computational 
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Figure 6: Form of the self-consistently determined voltage drop SUcs{x,y) averaged over 
the transversal y-direction for four different values of Ep . The results shown were obtained 
for the system investigated in Fig. [5] at bias Uq = 0.24 meV. For clarity, the curves are 
vertically offset by a value of —O.BUq each. For comparison, the dashed line shows the 
linear voltage drop considered in Fig. [5] 



feasibility. For a fixed bias Uq, it is given by 



In Fig. [7]we present the spin ratchet conductance (/^) /Uq as a function of the 
scaled driving amplitude IIq for three representative Fermi energies. In view 
of Eq. ( JT5l) . we expect the spin ratchet mechanism to be enhanced for higher 
driving amplitudes Uq, since the difference between P{+Uq) and P{—Uq) de- 
pends on the slope of the voltage drop. Indeed, in Fig. [7| we find a linear 
increase of the spin ratchet conductance for small Uq, which can be under- 
stood by expanding Eq. f|T5|) in this limit. 

We model impurity scattering by introducing Anderson disorder in the tight- 
binding version of the Hamiltonian ([3]) . To this end we add a random poten- 




for X < -NbLb/2 
for |x| < NbLb/2 
for X > NbLb/2 



(26) 
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Figure 7: Spin ratchet conductance {I^)/Uo as a function of the driving amphtude Uq 
for three different scaled Fermi energies, E-p = 0.24 (black lines), Ep = 0.33 (red lines), 
Ep — 0.44 (blue lines). Results are shown for a clean system (solid lines) and for systems 
with Anderson disorder of strength Udis = 0.5 (dashed lines) and Udis = ■\/2/4 (dotted 
lines). Parameters: Nb = 5, Lb = 10a, ao = 0-15, W ^ Wq = 15a, Ub = 0.22. 



tial from the box distribution [— f/dis/2; f/dis/2] to the on-site energy of each 
site in the central region of the wire. In Fig. [7] we compare the case of the 
clean quantum wire (solid lines) with results for the spin ratchet conductance 
for a disordered quantum wire (dashed and dotted lines). Although {I^)/Uo 
is reduced compared to the disorder-free system, the surviving spin ratchet 
conductance still possesses a reasonable magnitude. 

In order to get an idea of the corresponding mean free path in those calcula- 
tions, we now exemplary set the lattice spacing to a = lOnm and choose InAs 
with m* = 0.024 mo as the 2DEG material. Then the width of the ratchet 
device is W = Wq = 150 nm, the length of a single barrier is Lb = 100 nm 
and hence the overall length of the ratchet set-up is of order 1 micron. The 
Rashba SO strength is ao ~ 4.76 ■ 10~^^ eVm, a value well in reach of present 



day experiments |49|. The mean free path in this system can then be ap- 
proximated as / = 48aA/E^/f7|jg [sof. Then the disorder strengths chosen in 
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Fig. [TJ f/dis = 1/2 and correspond to an elastic mean free path of the 

order of one and two microns, respectively, for all Fermi energies considered. 
We hence can conclude from Fig. [7] that the spin ratchet effect prevails if I 
is of the order of or larger than the ratchet system size, i.e. even in the dis- 
ordered but non-diffusive limit. Since mean free paths I ^ l/xm are possible 



in clean InAs quantum wells |5ll], in realistic experimental situations spin 



ratchet output signals of reasonable magnitude should be observable. 
5. Conclusions 

In this paper, one the one hand, we have considered the spin- dependent 
transmission through single and multiple barriers in stripes build from two- 
dimensional electron gases with Rashba SO interaction. We have shown that 
for coherent transport across a smooth barrier with SO-mixed transverse 
channels, the spin transmission, i.e. difference between transmitted particles 
of opposite spin direction, increases linearly with channel number n. This 
implies an in-plane polarization of transmitted charge carriers, which is in- 
dependent of n or the wire width, respectively, and thereby distinctly larger 
than usual mesoscopic spin transmissions, such as e.g. from conductance fluc- 
tuations IJ, |52| or in the context of the mesoscopic spin Hall effect 53| which 
are typically of order 1. 

On the other hand, we have addressed the SO-mediated spin ratchet 
mechanism and thereby extended previous work, which gave a proof of prin- 



ciple for SO ratchets pj|, to realistic mesoscopic systems by accounting for 
certain nonequilibrium and disorder effects. In particular, we presented a 
self-consistent treatment of the ratchet transport in the strongly nonlinear, 
though coherent regime, showing that linear voltage-drop models often rep- 
resent good approximations to the real charge rearrangement but may break 
down in parameter regimes dominated by resonant charge and spin transfer. 
We further investigated the effects of elastic impurity scattering due to static 
disorder. We found that, as expected, the ratchet spin current decreases with 
increasing scattering rate, but a finite fraction of the clean spin current pre- 
vails (with the same sign) if the elastic mean free path is of the order of or 
larger than the system size. 

We close with a few further remarks: For systems with bulk inversion 



asymmetry, such as GaAs, Dresselhaus SO interaction [5J] has additionally 



to be considered. Calculations 36 for a ballistic SO ratchet with combined 



Rashba- and Dresselaus SO coupling showed that the overall picture is not 
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altered and furthermore demonstrate that the ratchet spin current direction 
can be changed upon tuning the relative strength of the two SO coupling 
mechanisms. In particular one finds spin current reversals for equal Rashba- 
and Dresselaus SO interaction where the SO effects cancel. This feature 
allows for inverting the polarization direction of the output current by simple 
electrical means, namely by tuning the Rashba SO strength, e.g. through a 
gate voltage. 

By now we have reached a fairly complete picture of the SO-based spin 
ratchet mechanism in both limits, the quantum coherent and the strongly dis- 
sipative regime. A further extension of previous work and a future challenge 
would consist in bridging these two separate limits treated up to now. 

In the coherent regime we so far considered the case of adiabatic driving 
that is relevant, e.g., for ratchet experiments employing nanostructures with 
an external ac bias voltage. If the system is driven through external radi- 
ation, the timescales for driving and for electron or spin dynamics can be 
comparable, which may give rise to further interesting spin ratchet phenom- 
ena. This ac regime requires an approach beyond the adiabatic limit, e.g. a 
Floquet treatment. 

The concept to extend the particle ratchet mechanism to other quantities 
such as the spin degree of freedom may be further generalized. Along this line 
one may think of extending this concept to other spin 1/2-type quantities, 
for instance the pseudo- or valley-spin degree of freedoms of charge carriers 
in graphene. More generally one may think of devising ratchet mechanisms 
to spatially separate objects according to their further internal degrees of 
freedom, for instance atoms with internal two- or multi-level dynamics. 
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